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■ The cross section and analyzing power of the pp pmr^ reaction near 

o : 

■ threshold are estimated in terms of data obtained from the pp dn^ and 
^ ■ 

r~| . pp ppn^ reactions. A simple final state interaction theory is developed 

. which depends weakly upon the form of the pion-production operator and 

^ : 

Ch ■ includes some Coulomb corrections. Within the uncertainties of the model and 

> '■ 

• "-^ . the input data, the approach reproduces well the measured energy dependence 

X: 

. of the total cross section and the proton analyzing power at a fixed pion cm. 

angle of 90°, from threshold to Tp = 330 MeV. The variation of the differential 

cross section with pion angle is also very encouraging. 
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I. INTRODUCTION 



The pp diT~^ reaction has been measured for many years [Q], either for its own intrinsic 
physics interest or as a cahbration for other reactions. However, because of the importance of 
estabhshing the precise beam energy, accurate experiments near threshold could be carried 
out much easier with the advent of cooled proton beams 0,^. The region within a few 
MeV of threshold is very interesting because it shows the interplay of the dominant p-wave 
production, induced by the A resonance, and the small s-wave rescattering terms. 

Since the phase-space for producing a three-body final state varies like Q'^, where Q is the 
cm. kinetic energy in the exit channel, the determination of the beam energy is even more 
crucial in such experiments near threshold. The accurate measurements of the pp — >■ pp7r° 
total cross section at the Indiana and CELSIUS storage rings down to Q ~ 1 MeV 
provoked great theoretical interest because the comparatively large cross section revealed 
the influence of unexpected short-range physics p. 

The total and differential cross sections for charged pion production in the allied pp 
pmt^ reaction have also been measured very close to threshold and these results have 
recently been complemented by the determination of the proton analyzing power . At the 
same value of Q, charged pion production is typically a factor of five larger than neutral 
and the interesting practical question to ask is to what extent the extra information will 
constrain theory further. 

The general features of the first results on the pp pmi^ total cross section near 
threshold [0] can in fact be understood within a nucleon-nucleon final state interaction 
model using measurements of the pp dir^ and pp ppn^ reactions as input |[10|. This 
approach uses an approximation for the np spin-triplet scattering wave function in terms of 
the deuteron bound-state wave function which does not depend upon the details of the np 
potential [|lT],|12|. This result, which is exact when the np scattering energy is extrapolated 
to —e^, where is the deuteron binding energy, is generally a very good approximation for 
small energies and distances [|T^ . Taking the pion production operator to be of short range, 
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the formalism allows one to estimate the cross section for the pp prni^ reaction where the 
final np pair are in the spin-triplet S'-wave. The much smaller contribution from spin-singlet 
final states may be deduced directly from the pp —>■ ppir^ data using isospin invariance. 

The agreement with the measured total cross sections is generally very encouraging 
except for underestimating the point nearest threshold [|T^ , which is the hardest to measure 
experimentally. In the definitive data set of Ref. p|, the value at this point is somewhat 
reduced and a more careful evaluation of Coulomb effects given in this present work increases 
the near-threshold prediction. These two changes reduce the discrepancy significantly. 

Away from threshold our final state interaction approach is capable of describing the 
angular distribution of the TRIUMF pp — > pmc^ data at 400 and 450 MeV [1^ provided 



the excitation energy in the np final state is below about 10 MeV, where S'-waves can be 



expected to dominate [|T5[. Differential cross sections in the pion cm. angle are now also 
available near threshold P], and we here extend our formalism to describe such data. A 
major uncertainty in implementing this program is the restricted differential cross section 
information for the pp — » ppir^ reaction corresponding to the total cross sections reported 
in Ref. ||^. Assuming this to be isotropic, we predict the pp prni^ angular distribution 
to be a little steeper than that observed experimentally though it should be noted that 
there are also uncertainties in the pp dir^ input . 

The proton analyzing power in low energy pp dn^ is largest at 90° and the energy 
dependence of this quantity for the three-body pp prni^ final state has just been reported 
by an lUCF group 0. The extension of our FSI method to treat analyzing powers is 
straightforward, but its application is hindered by the dearth of Ay data in the pp — > ppr:^ 
spin-singlet final state. Assuming this to be small, the broad features of the 90° energy 
dependence are reproduced. 

Calculations within our final state interaction model are not definitive in that they do 
not use all of the energy dependence of the S-state neutron-proton scattering wave functions. 
This means however that they do not require detailed considerations of the range dependence 
of the pion production operator. Higher partial waves in the neutron-proton final state are 



simply ignored and this will become a bigger defect with increasing energy. Nevertheless 
we claim that any microscopic calculation which reproduces the measured pp dn^ and 
pp ppiT^ observables should, near threshold, give results qualitatively similar to the ones 
presented here. The tolerable agreement with the pp — > pmi^ data therefore suggests 
that the extraction of quantitatively new information on the basic pion production reaction 
mechanism near threshold requires accurate measurements combined with much more refined 
theoretical calculations. 

In Sec. 2 we discuss the relation of the A^A^ scattering and bound state wave functions in 
the spin-triplet S-state and illustrate it with the results obtained from the Paris potential 
T3| . The energy dependence of the spin-singlet wave functions is investigated with and 



without Coulomb effects and it is shown that the major difference between the two can be 
accounted for by merely shifting the position of the ^Sq pole on the second sheet. 

Though at low energies it is reasonable to consider only S-wave nucleon-nucleon systems, 
pionic p-waves enter quickly because of the smallness of the s-wave terms. In Sec. 3 we 
develop the formalism to include the effects of pionic s- and p-waves, together with the 
approximate wave functions, in the three-body phase space to estimate the measured pp — >■ 
pmx^ observables. Coulomb corrections are included in the s-wave in the same approximate 
manner as that used to treat the low energy pp dix^ reaction. 

As previously stressed, the agreement with the experimental data presented in Sec. 4 is 
generally good, though not perfect. However the differences are well within the uncertainties 
of the model and the input and output experiments. Conclusions and suggestions for further 
work are to be found in Sec. 5. 

II. SCATTERING AND BOUND-STATE WAVE FUNCTIONS 

An S'-wave bound state reduced wave function behaves asymptotically as 

u^{r) ~ e-"'- (2.1) 
at large radii r and is normalised such that 



dr [u^{r)Y = 1 . (2.2) 

The analogous scattering wave function with real boundary conditions is normahsed by 
its asymptotic behaviour 

v{k,r) r^\si\i{kr + 5ik)) , (2.3) 
k 

where 5{k) is the phase shii^ at wave number k. 

Despite these very different normahsations, we have shown that for finite range potentials 



the two types of wave functions are related quantitatively by the theorem [T^ 



lim \j2a{a^ + k'^)v{k,r)\ = -u„(r) . (2.4) 

fe— >ja ' J 

This theorem is valid at the bound state pole, independent of the details of the potential, 
but it is also a robust extrapolation in the sense that, in the case of a weakly bound state, 
at short distances r 

v{k,r) ^ , ^ =Ua{r) , (2.5) 

^2a{a^ + k^) 

provided that k is small on the scale of the potential range and strength. Of course the ap- 
proximation breaks down at large distances because the bound and scattering wave functions 
must be orthogonal when integrated over all space. 

In the case of an almost bound virtual state, where the pole is on the second sheet, 
there is no bound state wave function to fix the radial dependence. Nevertheless the same 
techniques show that the energy dependence of the scattering wave function is given by a 
similar factor to that of Eq. ( p.5|) |[TO|.T^: 



^^^'^)^-^Wff^W- (2.6) 

These approximations can be tested explicitly in the cases of the Yamaguchi, square 
well, Bargmann or Hulthen potentials, which can be resolved exactly. However several 
complications arise when looking at more realistic potentials which describe the nucleon- 
nucleon system. The extrapolation theorem of Eq. (|2.5|) has only been proved for single 



channel scattering in a local potential [T^, whereas in the spin-triplet case there is coupling 



between the 5* and D-waves. 

In Ref. [|I0| we showed the variation of the deuteron and scattering wave functions with 



within 5% for r < 1.7 fm for neutron-proton cm. energies Enp < 10 MeV. To see this more 
quantitatively, we plot in fig. 1 values of the spin-triplet function 



with at = 0.232 fm~^, in terms of [|l^. The radius is fixed at r = 1.05 fm, which is close 
to the peak in the Paris wave function. 
The smooth fit 



lies 0.5% below the value of the deuteron wave function, but that is close to the limit in the 
precision of the extrapolation. This deviation could of course be due to effects of the D-state, 
though the Paris potential also includes some velocity dependence which may lie outside the 
domain of validity of the theorem. While almost all the energy dependence is given by the 
square-root factor in Eq. (|2.6|) , Ztik"^) is a steadily increasing function and this is typical 
for potentials with just one lightly bound state. Nevertheless, even at E^p = 15 MeV it lies 
less than 8% above the deuteron point. It should be noted that the slope of this function 
depends upon the value of r and that at r ^ 1.7 fm the function is almost fiat. 

Though there are no channel-coupling problems in the spin-singlet case, the approxima- 
tion given by Eq. ( |2.6|) must break down at very small k'^ in the proton-proton case due to 
the essential singularity at A; = caused by the Coulomb repulsion. This is not a problem 
provided -Eatat > 1 MeV, as can be seen from the values of the singlet function 



r and, despite the neglect of coupled channel effect; 



s, Eq. ( p^ ) seems to be broadly valid to 




(2.7) 




(2.8) 




(2.9) 



which are also shown in fig. 1. 
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The results may be parametrised by 

Z, = 0.646 + 0.343fc2 - 0.353fc^ + 0.157fc^ , with = -0.100 fm'^ (without Coulomb), 

= 0.698 + 0.131/^;2- 0.037/c^ with a, = -0.053 fm"^ (with Coulomb). (2.10) 

III. AMPLITUDES AND OBSERVABLES 

Many of the elements of the FSI formalism are given in Ref. |T^, but for clarity some 



are repeated here. We use a normalisation such that the two-body spin-averaged amplitude 
and cm. differential cross section are related by 

where mp, rrin, ma, and are the masses of the proton, neutron, deuteron, and pion, and 
and Pp are the momenta of the final pion and the initial proton respectively. 
The corresponding three-body cross section is given by 

-{pp ^ pn7r+) = — — -— X \M{pp ^ (pn)fc7r+)|2 , (3.2) 



d^l dk 32'K^{'mp + rrin + rriT^Ypp 

where the matrix element squared is assumed to be averaged over the angles of the np 
relative momentum k. In practice we shall only estimate the contributions from S'-wave np 
pairs, where the results are in any case isotropic in k. 

Reduced mass factors for the three-body reaction are defined by 

111 111 , , 

- = — + — - — ; — = — + — , 3.3 

/i TTLp + m„ rripn nip rrin 

whereas for bound state production the nip + m„ combination in Eq. ( |3.3| ) is replaced by 
the deuteron mass. 

The momenta in the exit channel are linked to the value of the excitation energy Q 
through 



The variable 77, conventionally used to describe low energy pion production, is the maximum 
pion momentum in pion mass units and is given by 

r] = = . (3.5 

The cross section for pion production summed over all pn excitation energies in the 
three-body channel is 



do 1 r^'irripriQ 

— {pp^pn7c^) = —— ■ ■ — / \M{pp ^ {pn)k7i^)\ p^k dk . (3.6) 

d\l S2TT'^[mp + TUn + m-^j^pp Jo 

Since the pion production operator is expected to be dominated by short-range physics, 
the relation of Eq. (|2.4|) between the bound and scattering wave functions suggests a similar 



relation for the spin-triplet matrix elements, viz 

\M{pp ^ {pn)t7i^)\ ^ ^ ^ —{pp dTi+) . 3.7 

2at[k^ + at) pT,at[k^ + af) d\l 

The prediction for the spin-triplet part of the unpolarised cross section then follows imme- 
diately through the use of Eq. ( |3.6| ) . 

The formalism for the proton analyzing power Ay in pp pn7[~^ can be developed 

identically, the only changes being Ay factors on the right hand side of Eq. ( |3.71 ) and the 
left hand side of Eq. ( |3.6| ) . 

It is of great practical interest that near threshold the integration in Eq. ( |3.6| ) may be 
performed analytically. Thus for a pp dn^ observable which varies as 

0^;\V)^V\ (3.8) 
the corresponding observable for the pp — >■ {pn)tn~^ reaction is given by 

(^i"J(^) = ^^"H^C)xOl")(77), (3.9) 

where 

C = ^- 3.10 
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A description of the low energy pp prni^ observables only requires these functions for 
n = 1, 2, 3, and explicitly in these cases 



P^^\x) = " (3.11a) 

^ ' 4(l + v^rT^)2 ^ ' 

p(^\x) = -l-x + -- + arctan(x)l , (3.11b) 

TT 13 X \ X^ } J 

P''\^) = ... . I . \ 1 + ... . I . ,,,, \ ■ (3.11c) 



[i + VT+^y\ 2(1 + va 



X 



It is customary in the analysis of pp dvr^ data to assume that, due to the Coulomb 
repulsion in the final state, the cross section very near threshold is suppressed by a Gamow 
factor of 

2q/ri 1 
^^^''^ = exp(2g/r^) - 1 ^ TT^ ' ^^'^^^ 

where 

q = (3.13) 

This approach, which considers all the charge in the deuteron to be concentrated at the 
center of the nucleus, is sufficient to explain the major differences between the np — dn^ 
||T6| and pp —>■ dn^ ||T|] total cross sections, and the simplification made in Eq. ( |3.12[ ) is 
indistinguishable from the full Gamow factor under the conditions of such data. 

An analogous suppression will be suffered by the tt"*" in the pp —* pmx^ case, but this 
is only likely to be of practical importance for s-wave pions because the centrifugal barrier 
keeps the p- waves small in the region of low rj. Taking the charge, as in the deuteron 
case, to be situated at the center of mass of the proton-neutron system, the integral for the 
Coulomb-modified P^^-* {x) of Eq. ( |3.11|) can still be performed in closed form, giving 



9 



3 r /i_^/r^72 



X 10 n ^ „ / i — V i — C 



+ ^(1 + x\l - c^m + x\l + 2c^)) - ^ (1 + x^)^/^ 

-4(1 + x^n - c^)) + 4 (1 + x^) arctan(x)l , (3.14) 

X"^ X'^ J 

where c = q/rj. 

In pp — ^ ppn^ the s-wave pion production is hkely to be dominant to much higher values 
of T] than for charged pion production induced by the suppression of the A contributions 
through the selection rules for the final S'-wave pp states. Though there is no bound state 
production to set the normalisation, the energy dependence of the total cross section should 
therefore be of the form r]P^^\x). In the absence of Coulomb effects, the cross sections for 
7r° and production are related by 

a{pp {pn}/=i7r^) = a{pp ppn^) . (3.15) 

Coulomb corrections must however be applied before this is used to estimate the spin- 
singlet contribution to pp pn-K^ . On the right hand side the singlet pole in the wave 
function is shifted as in Eq. ( |2.1CI|) , whereas on the left the Coulomb modified P^\x,c) is 
used instead of P^^'^x). 

IV. COMPARISON WITH EXPERIMENT 

In the low energy region it is a good approximation to take the pp d7c~^ differential 
cross section and analyzing power to be of the form 

^ = A + Bcos^e, (4.1a) 
Ay^ = Csm9, (4.1b) 

where 6 is the cm. angle of the pion, and these dependences will propagate through to the 
spin-triplet part of the pp pmT~^ reaction. 
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The SP96 phase shift solution from the Virginia database yields the following 
Coulomb-corrected predictions for observables. At 6* = 0°, 

d(j 



(13.6 7/ + 164.0 r?'^) /ib/sr, (4.2) 



whereas at 6' = 90° 



(13.6 T] + 16.2 T]^) ijh/sT , (4.3a) 



dn* 



J^ = -18.1r/Vb/sr, (4.3b) 

though it should be noted that the s-wave strength reported here is 7% less than some direct 
measurements |0,|16 . 



Spin-singlet final states may be estimated using 

a{pp pjM^) ^ 4nD r] P^^^ {('r]) (4.4) 

where, as discussed in Sec. 2, after including the pp Coulomb repulsion, (' = 12.6. The best 
fit to the total cross section data of Ref. is achieved with D ^ 0.86 /ib, and the resulting 
curve shown in fig. 2 is very similar in shape to microscopic calculations which include the 
pp final state interaction correctly UTTf . 

Though the authors of Ref. did not report values for the pp —>■ ppn^ differential cross 
sections, their data do suggest some influence from higher partial waves at 325 MeV. There 
is however evidence that at 320 MeV (?7 = 0.46) the cross section is fairly isotropic and that 
the proton analyzing power is consistent with zero [|T^ . Introducing Coulomb corrections in 
both the final ppn'^ and ppn^ systems, these data then predict that 

(pp ^ {pn}sn^) ^ D'r^P^'\C\c) , (4.5) 



dn* 

where D' ^ 0.42 /ib/sr and (" = 25.4. This is to be taken in conjunction with the spin-triplet 
predictions for the differential cross section 

■^{pp ^ {pn}t7r^) = 13.6 r]P^}\C'v,c) + (16 + 148 cos^ ^) r/^ P(^)(C"t/) /ib/sr, (4.6) 
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and analyzing power 

Ay -^{PP {pn}t''T+) = -18r]'^smeP^^\C"r]) fih/si , (4.7) 

where pion Coulomb corrections have only been made in the s-wave. 

The momentum dependence of the pp — > pnir^ total cross section is shown in fig. 2 
both with and without pionic Coulomb corrections. There is in fact very little difference 
between either using the P(^^''(x,c) of Eq. (|3.14|) or multiplying the uncorrected P^^\x) of 



Eq. (|3.11| ) by the Gamow factor of Eq. (p.l2|) . The full curve lies somewhat below the 



lUCF measurements § but, apart from the lowest point, the deviations are of the same 
order as the Coulomb differences. It should be noted that there are typically 10% overall 
normalisation uncertainties in both the input and output data, as well as uncertainties in 
determining the value of rj close to threshold p[. The total cross sections from the extended 
data set of Ref. are consistent with the earlier measurements, though they are subject to 
very similar systematic uncertainties. 

Due to the differences in the pole positions, the ratio of the s-wave production of pions 
with spin-triplet and spin-singlet final states is estimated to be about 1.7 at low rj but 7.8 
at high 7]. 

The angular dependence of the predicted pp prnr'^ differential cross section at = 0.21 
and rj = 0.42 is compared in fig. 3 with the lUCF data ||^. The theory is about 20% too low 
in normalization and the slopes, evaluated assuming the singlet cross section to be isotropic, 
a little too strong. It is clear from the larger data set of Ref. that the binning effect due 
to the experimental averaging over angular domains is only of marginal importance. 

In Fig. 4 is shown the variation of Ay at QO'' as a function of rj for pp — > dn^ and 
pp pmT~^, assuming the singlet cross section to be isotropic with zero analyzing power. 
The curve lies a little below the experimental points of Ref. [^j, where the only major 
systematic uncertainty is the acceptance in polar and azimuthal angles. 
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V. CONCLUSIONS 



We have shown within the framework of a simple final state interaction theory that most 
of the low energy data on the total and differential cross sections and proton analyzing power 
for the pp pmi^ reaction can be understood semi-quantitatively in terms of equivalent 
information deduced from pp —>■ ppn^ and pp —>■ dn^ measurements. We would argue 
that any more microscopic approach with finite range pion-production operators, which fit 
the latter data, would give pp — *• pmr~^ predictions broadly similar to the ones presented 
here, though with some deviations at higher values of rj due to the neglect of final np P- 
wave contributions. Such contributions might be responsible for some of the discrepancies 
observed in the angular distributions and analyzing power. To take full advantage of the 
accurate lUCF data would require a sophisticated microscopic model to be accurate to the 
10% level. It could then be helpful if such a model were consistently applied in the evaluation 
of the experimental acceptances rather than relying on the Watson FSl approach P]. 

Qualitatively new information, against which to test our approach, would be provided by 
experiments which introduce interferences between triplet and singlet final np states, as for 
example in the {p, p') spin-transfer. These are however much harder to perform in practice. 

The deviations between our predictions and the lUCF data [0-3 may be partially ex- 
perimental in origin. It should be noted that there are uncertainties of the order of 10% 
in some of the pp dTT~^ observables used as input, and this is worse for the differential 
quantities in pp ppn'^. However it is also apparent from fig. 1 that our formalism will 
tend to underestimate slightly the triplet contribution to the pp prnr^ cross section be- 
cause Zf is an increasing function of k"^. The rate of rise depends upon the nucleon-nucleon 
separation and so it is impossible to quantify the effect without having a detailed model for 
the pion-production operator. This is therefore left for more microscopic formulations. 
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Figure 1: The NN triplet scattering function Zt defined by Eq. ( p.7|) evaluated at r = 1.05 fm 
for discrete values of k"^ using the Paris potential (open circles). The smooth extrapolation 
of Eq. ( p.8|) lies slightly above the deuteron point (closed circle) . The points for the corre- 
sponding singlet functions Z^, defined by Eq. ( ^.91 ), are shown without Coulomb force (open 
star) and with (closed star) and interpolated using Eq. (|2.10|) . The principal effect of the 



Coulomb force is to change the position of the antibound state pole at k"^ = —a^. Residual 
effects due to the Coulomb repulsion may be seen below Epp = 1 MeV. 
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X! 




Figure 2: The experimental pp — ppn^ total cross sections of Ref. 0] (open circles) and Ref. 
(squares) are to be compared to the broken curve, which represents the fit on the basis 
of Eq. 
The pp 



This curve is very similar in shape to the microscopic calculations of Refs. [|17 
data of Ref. 



pnvr ' aata ot Jrtet. lajj, which are very similar to the extended data set of Ref. 0, 
are subject to a 13% overall uncertainty. The solid and chain curves show the corresponding 
predictions with and without the pionic Coulomb corrections. The Coulomb corrections 
given by Eq. (3.14) differ little in practice from those obtained by merely multiplying by the 
Gamow factor of Eq. (3.12). 
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Figure 3: Differential cross sections for pp — »• pnir^ at r] = 0.21 (squares) and rj = 0.42 
(circles) taken from Ref. [§] and compared with the predictions of the model assuming the 
singlet cross section to be isotropic. The data and predictions at the lower energy are both 
multiplied by a factor of five to present them on a similar scale. 
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Figure 4: Variation of the proton analyzing power in pp — > prirr'^ at 6 = 90" as a function of 
rj, the experimental points being taken from Ref. 0. The broken curve is the prediction of 
the phase shift solution SP96 of Ref. [||] for pp dn^, whereas the solid curve follows from 
Eq. ( [4.7| ) assuming the pp — > ppvr" cross section to be isotropic with neghgible analyzing 
power. 
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